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Abstract. I review recent developments in finite temperature lattice QCD which are 
useful for the study of heavy-ion collisions. I pay particular attention to studies of the 
equation of state and the light they throw on conformal symmetry and the large N c 
limit, and to the structure of the phase diagram for Nf = 2 + 1. 
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1. Introduction 

Lattice gauge theory at finite temperature and density has been an extremely active 
field in the last year. This makes it hard to review it in its entirety within the scope of 
this article. I have therefore chosen to review two of the major points of contact with 
heavy-ion physics, namely the equation of state and the phase diagram, in the later 
sections. In this section I provide pointers to the literature on the many developments 
I will not discuss further. 

• The finite-temperature phase transition is now definitely established to be a cross 
over: this is verified by the Budapest- Wuppertal (BW) group [TJ[2], and confirmed 
by the BBRC collaboration [3j 0]. Further confirmation has come from the Hot- 
QCD collaboration. The cross over temperature is temporarily in dispute. The 
Budapest- Wuppertal group [2] finds a cross-over temperature substantially smaller 
than that obtained by BBRC [4j using the same thermometer. An old global 
analysis in 2001 [5] gave T c ~ 175 MeV with 20 MeV uncertainty from scale 
setting. BBCR and Hot-QCD prefer the upper end, BW prefer the lower end. This 
disagreement, while serious, therefore does not impact the experimental search. 

• Correlation function measurements have led to new developments. Deconfinement 
occurs at the chiral cross over point: this has been shown using linkages 
between quantum numbers; for example, the linkage between baryon number and 
strangeness becomes exactly that expected from quarks at T c [6]. There is steady 
and slow advance in the measurement of transport coefficients on the lattice [7\. 
Renormalized Polyakov loop measurements in various representations give strong 
evidence for Casimir scaling to all orders [8]. The octet loop, in particular, does 
not see the phase transition. The non-melting of J/ip and the melting of Xc soon 
above T c is now verified in many different computations [9]. 

• There has been significant advance in defining chiral fermions at finite chemical 
potential [10J. There are advances in the understanding of isospin chemical potential 
and imaginary chemical potential [IT] . There are new results for thermodynamics 
using Wilson quarks [12]. Localization of staggered Dirac eigenvectors is seen to 
set in abruptly at T c , although this could be a finite volume artifact [T3] . 

2. The equation of state 

The equation of state is rather well-known for QCD without quarks. Recently there have 
been three separate efforts to extract it for Nf = 2 + 1, i.e., two light degenerate flavours 
(corresponding to u and d quarks) and one heavier flavour. The BW collaboration uses 
staggered quarks and sets the renormalized quark masses to their physical values by 
tuning to realistic values of and mx- The BBRC collaboration uses P4 quarks with 
light quark masses heavier than physical but a physical strange quark masses by tuning 
to a heavier pion (m n ~ 220 MeV) and a realistic m^. The HotQCD collaboration is 
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Figure 1. Quenched data from |14| . Nf = 2 + 1 data from [2JH]. On the left, is a plot 
of P/T 4 against E/T 4 . On the right is a plot of P/P SB and E/E SB - The quenched 
data is scaled by the Nf = continuum SB values, Nf = 2+1 data is scaled by the 
Nf = 3 continuum SB values. 

investigating the equation of state with both P4 and Asqtad quarks at renormalized 
quark masses equal to that used by the BBRC collaboration. 

The equation of state of state in the form P(E) is an important input into 
hydrodynamics, and therefore very important for predictions of various signals in heavy- 
ion collisions. In conformal field theories, including free field theory, one necessarily has 
P = E/3, so that c 2 s = 1/3. This limit has excited interest recently, since many toy 
models for QCD can be solved using AdS/CFT techniques which demand conformal 
symmetry along with the large N c approximation. 

The lattice data is displayed in the two composite plots in Figure HJ The diagonal 
line denotes conformal equations of state: points corresponding to the ideal gases for 
Nf = 0, 2 and 3 are marked on this line. Note that the data on the quenched theory 
is extrapolated to the continuum limit, whereas those for Nf = 2 + 1 are for a finite 
lattice cutoff a = 1/4T and a = 1/6T. There is apparent disagreement between the 
computations with staggered and P4 quarks: they are not expected to agree at finite 
lattice spacing, only in the continuum limit must they give the same result. In this 
connection note that the results of the HotQCD collaboration at a = 1/8T smoothens 
the behaviour of the P4 quarks, reducing the maximum excursion from the conformal line 
by about 20%, while leaving data at 3T c /3 and higher almost unchanged. In summary, 
there is clear evidence for strong violation of conformal symmetry in the region of the 
crossover from hadron to quark phases. 

Importantly, the deviations are not restricted to a narrow region around T c , but 
extend to T = 2T C or slightly higher. This will have two consequences in applications 
to heavy-ion collisions — 

(i) Hydrodynamics of non-ideal fluids with conformal equations of state (P = E/3) 
allow only for shear viscosity. If conformal symmetry is broken, the fluid can have 
bulk viscosity as well. Closely related to this is the fact that c 2 s < 1/3. 

(ii) Unless initial temperatures rise well above 2T C , heavy-ion collisions may never be, 
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Figure 2. Strong deviations from conformal symmetry are observed for N c = 4 and 
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a priori, approximated as systems with vanishing bulk viscosity. 

A recent study of SU(A r c ) pure gauge theories shows that this behaviour persists 
for N c > 3 [16J. The breaking of conformal symmetry is often exhibited as a plot of 
E — 3P against T (see Figure [2]). For SU(3) gauge theory with or without quarks, peaks 
in this quantity are seen close to T c . For pure gauge theories, there is a first order phase 
transition, and the maximum value of E — 3P is at least as large as the latent heat 
density, and hence expected to scale at N%. Since the Stephan-Boltzmann values of E 
and P also scale as N c , the breaking of conformal symmetry for SU(iV c ) pure gauge 
theories is also expected to persist to large N c . 

In heavy-ion collisions a very important use of the equation of state is the estimate 
of typical energy densities in the plasma phase and in the hadronic phase of QCD. For 
reference I quote some estimates for the energy density from current lattice computations 
just below the lowest estimate of the crossover temperature, and just above the highest 
estimate. We can take as fiducial values E\ = £7(140 MeV) which is definitely in the 
hadron phase and E-i = £7(210 MeV) which is definitely in the plasma phase. The 
current estimates are 



70 ± 12 MeV/fm 3 (BBRC, 1/a =840 MeV), 

21 ± 24 MeV/fm 3 (BW, 1/a =840 MeV). 



In the high temperature phase we have, 



3.11 ± 0.07 GeV/fm 3 (BBRC, 1/a =1260 MeV), 

2.28 ± 0.08 GeV/fm 3 (BW, 1/a =1260 MeV). U> 

While noting the statistically significant differences between the results of the two 
computations, note also that the renormalized quark masses differ, and that both 
estimates are made at finite lattice spacing with quark formulations which differ at 
finite cutoff. For comparison, note that in the continuum limit of quenched QCD one 
finds E(T C = 285 MeV) = 3.4 ± 0.5 Gev/fm 3 [TJ]. 
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Figure 3. The phase diagram on the left is put together from computations by [TH] 
and arguments from |19j . Projecting the phase diagram down to the T = plane, one 
obtains the flag diagram, which is the figure on the right. 



The success of a hydro dynamic description of heavy- ion collisions would eventually 
be gauged from its ability to reliably extract P(E) (and other material properties 
such as transport coefficients) from data, thus permitting experimental tests of lattice 
predictions. Any claims of disagreement between such extractions and lattice predictions 
are actually claims of the failure of QCD, and therefore should be treated with the 
strong, fair and rational skepticism that any claim of a failure of a well-established 
theory should receive. At this point of time it is easier to control lattice studies: one 
merely has to spend computer time in decreasing the lattice spacing, a procedure that is 
already under way. Control of the hydrodynamic description is more problematic, with 
initial conditions, the hydrodynamic equations, and the hadronization prescription each 
still in need of independent validation. 

3. The global phase diagram 

Phase diagrams are labeled by the thermodynamic intensive coordinates. For QCD 
these are T, Nf quark masses and Nf chemical potentials. Of these, experiments can 
tune (at best) 1 + Nf of these, since the quark masses are fixed conditions that we 
are faced with in reality. In practice, heavy- ion collisions have only a single control 
parameter — y/S. This is not sufficient to examine the 4 dimensional phase diagram 
of QCD: only enough to explore a single line through the phase diagram. By varying 
the ions one can smear this line a bit, but this still leaves scope for much thought and 
experimental ingenuity in exploring larger parts of the phase diagram in the laboratory. 

Each point in phase diagram is, almost always, a single pure phase. Exceptions 
are where two or more phases coexist; these are also called first order transitions. A 
continuity argument for lines (surfaces) of first order transition is called the Gibbs' 
phase rule. It follows from the structure of the solutions of the equation gA(T, m, rrii) = 
gs(T, fiijUii), where qa,b are the free energy density in the two phases A and B. The 
Gibbs' phase rule implies that in D dimensional phase diagram one has D-2 dimensional 
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Figure 4. The extended flag diagram of QCD. The planes /i = and m s — oo of this 
flag diagram are reasonably well explored, and have the structure shown in the first 
panel. There are two tricritical points, at each of which an Ising critical line joins an 
0(4) critical line. Since the two 0(4) critical lines bound an 0(4) critical surface, the 
two tri-critical points are joined together by a tri-critical line. This line has to lie on 
the m — plane, since the the 0(4) surface lies in that plane. In [T7] it was found 
that a part of the Ising critical surface near the Nf = 2 + 1 region bends "backwards" . 
Since a critical surface does not have "surface tension" , it can twist and bend. To look 
at the possible phase diagrams, examine the shape of the tricritical line. 



critical surfaces, D-3 dimensional tricritical surfaces, D-4 dimensional tetracritical 
surfaces etc.. This argument, and equivalent forms of it, strongly constrains the topology 
of phase diagrams [20J. 

When discussing high dimensional phase diagrams, one often discusses sections, i.e., 
parts of the phase diagram with some of the intensive variables set equal to fixed values. 
These are also phase diagrams, in the sense that each point corresponds to an unique 
thermodynamic phase. However, in QCD it has become the practice to use projections, 
which we call flag diagrams. A construction of a flag diagram is shown in Figure [3j 
Each point in a flag diagram shows the kind of phase transition that one obtains by 
varying the quantity which has been projected away. 

The QCD flag diagram can be extended to include finite baryon chemical potential 
fi. Note that we are discussing a 3-d flag diagram obtained from the 4-d section of the 
full phase diagram with m u = = m, and the two chemical potentials, /i s = \ij — 0. 
In this flag diagram the planes of m s = oo and \i = have been explored, and the 
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Figure 5. Three topologies for the tricritical line are shown. Below each is shown 
the phase diagram corresponding to two different generic quark masses. Note that the 
second case above is ruled out by universality arguments. 



topology of the phase diagrams are reasonably well understood. In each of these planes 
there is a tricritical point at which an Ising critical line joins on to an 0(4) critical line. 
One of these tricritical points is shown in Figure [3) Now, the two 0(4) critical lines are 
just the boundary of a single 0(4) critical surface. Since these two end in a tricritical 
point, there must be a tricritical line bounding the 0(4) critical surface, of which these 
two points are the ends (see Figure HJ). 

Along this tricritical line one must glue together an 0(4) critical surface and one 
or more Ising critical surfaces. Now, [17] found that a part of the Ising critical surface 
near the Nf = 2 + 1 region bends in the direction of smaller m s . With this evidence 
they claimed that the tricritical line must occur in two pieces, and the two known Ising 
critical lines must belong to two different Ising critical surfaces. This argument, while 
compelling, is not water-tight. The reason is that there is no physical principle that 
constrains the curvature of critical (and tricritical) surfaces: they are allowed to bend, 
twist and wander. Hence the "wrong" curvature seen in [17] can be accommodated 
into the flag diagram shown in Figure HI One must enumerate all topologies that the 
tricritical line may have. 

Apart from the simple topology shown in the last panel of Figure HI there are 
several more possibilities. In a 4-d phase diagram one non only has the possibility of 
a tricritical line, but also a tetracritical point where two such tricritical lines meet. If 
indeed the two Ising critical surfaces are physically distinct, then the possibility of a 
tetracritical point must be considered. Usually, one would discover a higher-order critical 
behaviour at a point of enhanced symmetry. Such points are along the 3 flavour line 
m s = m„ = ma- However, this is ruled out, since universality arguments rule out critical 
behaviour for SU(3) flavour. The remaining possibility is to draw the tetracritical point 
out to infinity. These three distinct possibilities for the topology of the tricritical line] 
are shown in Figure along with the three dimensional physical phase diagrams that 
they would give rise to at generic quark masses. Note that the second case is ruled 
out by universality. Which of these cases is true for QCD is something that has to be 
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decided by further lattice simulations. 



4. Summary 

Significant progress has been made in lattice computations in the last year. The new 
state of the art is m n ~ 140-220 MeV. For the first time all collaborations are using 
< m p /2. Conformal symmetry is strongly broken for T < 2T C , as shown directly 
by the equation of state. There is direct evidence from the lattice that conformal 
symmetry breaking persists for larger number of colours, and hence, through power 
counting arguments, to all N c . One major consequence of this observation is that 
hydrodynamic studies should include bulk viscous terms. I have argued that a QCD 
critical point could exist at small chemical potentials even in the Nf = 2 + 1 theory, 
and that it is still possible that Nf = 2 computations are a good guide to its location. 
Further lattice computations are called for. 
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